Metal-insulator transition in graphene induced by circularly polarized photons 
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Exact stationary solutions of the electron-photon Dirac equation are obtained to describe the 
strong interaction between massless Dirac fermions in graphene and circularly polarized photons. 
It follows from them that this interaction forms bound electron-photon states which should be 
considered as a kind of charged quasiparticles. The energy spectrum of the quasiparticles is of 
dielectric type and characterized by an energy gap between the valence and conductivity bands. 
Therefore the electron-photon interaction results in metal-insulator transition in graphene. The 
stationary energy gap, induced by photons, and concomitant effects can be observed for graphene 
exposed to a laser-generated circularly polarized electromagnetic wave. 

PACS numbers: 73.22.-f 



I. INTRODUCTION 

Since the discovery of graphene^ — a monolayer of 
carbon atoms — its unique physical characteristics have 
aroused enormous interest in the scientific community. 
Particularly, the influence of electromagnetic field on 
electronic properties of graphene is in focus of attention". 
However all theoretical models, elaborated to describe 
this influence, deal with strong electromagnetic field in 
the framework of classical electrodynamics. A theory, ca- 
pable of describing the strong interaction between elec- 
trons in graphene and quantized electromagnetic field, 
was unknown up to the present. The paper is aimed 
to fill partially this gap in the theory. It is surprising 
that the problem of interaction between electrons at the 
Fermi level of graphene and circularly polarized photons 
can be solved exactly in analytical form without the need 
of any approximation and numerical calculations. This 
is one of few exactly solvable problems in the theory of 
photon-matter interaction, that is not without interest 
to both condensed matter physics and quantum optics. 
Moreover, since electrons in graphene behave as mass- 
less relativistic fermionsii^, graphene can serve as a 
low-energy 'proving ground' for the relativistic physics 5 . 
Therefore the obtained solutions can also be useful to 
describe the interaction between relativistic particles and 
photons. Thus the solved problem lies at the intersection 
of different excited fields of modern physics, that is what 
caused to write this paper. 

II. DIRAC PROBLEM FOR 
ELECTRON-PHOTON SYSTEM IN GRAPHENE 

Generally, electron states in graphene near the Fermi 
energy are described by eight-component wave functions 
written in a basis corresponding to two crystal sublat- 
tices of graphene, two electron valleys, and two orienta- 
tions of electron spin^. In what follows intervalley scat- 
tering processes and spin effects will be beyond considera- 
tion, that reduces the number of necessary wave-function 
components to two. To describe the interaction between 



a plane monochromatic circularly polarized electromag- 
netic wave and electrons in a graphene, we shall be to use 
the Cartesian coordinate system (x,y,z), where the axis 
z is perpendicular to the graphene fixed at z = 0. Then 
near the point at the Fermi level, where the valence and 
conductivity bands of graphene touch each other (the 
Dirac point), the Hamiltonian of electrons interacting 
with the electromagnetic field can be written for a single 
valley and for a certain direction of electron spin in the 
form 2 

H e = V F & (p - -A) , (1) 

where vf is Fermi velocity, p is operator of electron mo- 
mentum in the graphene plane, e is electron charge, and 
A is vector potential of the electromagnetic field. As to 
the vector operator <r, its components are Pauli matri- 
ces written in the basis of two orthogonal electron states 
arisen from two crystal sublattices of graphene. Since 
the operator a z is diagonal, these two states can be 
denoted by mutually opposite orientations of the pseu- 
dospin along the z-axis, s z = ±1/2. Considering the 
problem within the standard quantum- field approach 6 , 
the classical field, A, should be replaced with the field 
operator, A. Assuming the electromagnetic wave to 
be clockwise-polarized and propagating along the z-axis, 
this operator can be written as 

A = ^2Trhc 2 /u V (e+a + e_a f ) , (2) 

where ujq is frequency of the electromagnetic wave, V is 
volume of space with the field, e± = (e x ±ie y )/y/2 are po- 
larization vectors, e XjV are unit vectors directed along the 
x, y-axes, a and are photon operators of annihilation 
and creation, respectively, written in the Schrodinger rep- 
resentation (the representation of occupation numbers) 6 . 
Then the complete Hamiltonian of the electron-photon 
system, including both the field energy, hujoa^a, and the 
electron Hamiltonian ([1]) rewritten in the quantum-field 
form, is given by 



% = huod^a + vf0"P — e\ —f- (o+a + , (3) 

V 
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where a± = (a x ± ia y )/2 are step- up and step-down op- 
erators for the z-projection of the pseudospin. 

The stationary solutions of the effective Dirac equation 
with the Hamiltonian §3§ have the form 'f(k) = e lkr f/'(k), 
where k is wave vector in the graphene plane, r is radius- 
vector, and ^(k) is eigenstate of the Hamiltonian 



uj V 



<j + a + (7_a 



(4) 



At the Dirac point (k = 0) the Hamiltonian ^ is 
formally similar to the Hamiltonian of exactly solvable 
Jaynes-Cummings modeli. As a consequence, the effec- 
tive Dirac equation with the Hamiltonian (f4]) for k = 
can also be solved exactly. To describe the electron- 
photon system at the Dirac point, let us use the nota- 
tion \s z ,N) which indicates that the electron is in one 
of two quantum states with the pseudospin projections 
s z = ±1/2 and the electromagnetic field is in quantum 
state with the photon occupation number A = 1,2,3,.... 
Then exact eigenstates of the Hamiltonian (£[]) for k = 0, 
■0(0) = y+^jv an d V'(O) = y_i,jv) can be written as 



\<P±l, n) 



±1/2, TV) 



where 



16(A + 1/2 ± l/2)(irv 2 F e 2 /huJoV) + w, 



(5) 



energies of the electron-photon states ([5]), e , i N and 



are given by 



e ±%,N = Nhu ± — T ~y~ 



(6) 



and subscript indices in Eqs. ([H])-© indicate genesis of 
these states: The state \<p±i at) turns into the state 
| ± 1/2, A) when the electron-photon interaction vanishes 
(i.e. for e = 0). The expressions ©-(HI can be eas- 
ily verified by direct substitution into the effective Dirac 
equation Hk^ii n = n l P±\ n with the Hamiltonian 
(|3]) for k = 0, keeping in mind the trivial relations^ 



a±\Tl/2,N) = | ±1/2, AO, &±\±1/2,N) = , 
a f | ±1/2, N) = vf + T | ± 1/2, A + 1) , 
a|±l/2,A) = VN\±1/2,N-1). 

Let the photon occupation number of the electromagnetic 
field is fixed at A = Ao by a field source. Then the 
electron-photon states ip_i No and (p + i No , originated 
from electron states degenerated at the Dirac point, are 
separated by the energy gap e g = e_i_ No — e + i Nq . In 
what follows we shall be to assume the electromagnetic 



wave to be classically strong, that corresponds to macro- 
scopically large photon occupation numbers (Nq S> 1). 
Then the gap can be written as 



(7) 



where Wq = 2vfgEq/ujq is energy of electron rotational 
motion induced by the circularly polarized wave, and 
Eq = y^AirNohuo/V is classical amplitude of electric field 
of the wave. 

Let us proceed to solving the effective Dirac equa- 
tion with the Hamiltonian (j4|) near the Dirac point for 
wave vectors k / 0. Since the obtained electron-photon 
states ([5]) are eigenstates of the same Hamiltonian (f4]) 
for k = 0, they form complete basis of the considered 
electron-photon system. Therefore true eigenstates of the 
Hamiltonian (0} for k ^ can be sought as an expansion 



0(k) 



E 
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Substituting the expansion (j8]) into the effective Dirac 
equation Hkip(k) = e(k)^>(k) with the Hamiltonian (j4}, 
we arrive at the system of recurrent algebraic equations 
for coefficients c ± i m (k). This system can be easily 
solved in the special case, when the effective parameter 
of electron- field interaction, a = Wo/haj , satisfies the 
condition a -C 1. Seeking the dispersion dependencies 
e(k) and c±i m (k) for any k in the principal order with 
respect to a, the system of recurrent equations can be 
reduced to the two equations 

(e - NqHwq - ^ ) c_i jy = hv F (k T . + ik„ 



(e - N a fuv a + -Z 



-+h,No 



hv F (k x 
%v F {k x 



,n , (9) 

,N ,(10) 



which can be solved exactly. Solving Eqs. (l^- flTUl) by 
conventional methods, we obtain two branches of the en- 
ergy spectrum, e(k), which can be written as 



4 o (k) - N hw ± y / (e a /2)2 + (^ F fc)2, (ii) 

where omitted k-terms are ~ o(a). As a result, eigen- 
states of the Hamiltonian ([3]) are given by 



.N 



(k)v>- 



,N B 



(12) 

where superscript indexes '±' denote solutions of 
Eqs. (l9|- (rrtH) corresponding to the two energy branches 
(ITT1) . As expected, for k = the obtained expressions 
(TTT1) - (fT2"|) coincide with the exact solutions ©-([5]). 



III. ELECTRON-PHOTON QUASIPARTICLES 
IN GRAPHENE AND THEIR PROPERTIES 

The full energy of the electron-photon system (ITT1) con- 
sists of the field energy, NqHujo, and the term 



Ae(k) = ±./(e 9 /2) 2 + (hv F k) 2 , 



(13) 
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which is arisen from electron-photon interaction. As a 
consequence, an electron interacting with circularly po- 
larized photons can be considered formally as a quasi- 
particle with the energy spectrum (|13p . From the view- 
point of classical electrodynamics, the circularly polar- 
ized electromagnetic wave rotates the electron 9 -. There- 
fore the quasiparticle is the mix of photon states and 
electron states corresponding to this rotational motion 
of the electron. Certainly, a center of trajectory of the 
rotating electron can move in the graphene plane, that 
corresponds to a translational motion of the quasipar- 
ticle as a whole. The translational motion is described 
in Eq. (fT3]) by the wave vector of quasiparticle, k. Since 
the quasiparticle is rotating electron dressed by circularly 
polarized photons, this wave vector differs from the wave 
vector of free electron in graphene. 

The energy spectrum of quasiparticles <p~3|) . schemat- 
ically pictured in Fig. la, defines optical and transport 
properties of graphene irradiated by circularly polarized 
light. The distinctive feature of the spectrum (|T3"|) is the 



Ae(k)' 

\ 


As (k)' 




> w 




\ 




> 
> 


/ j# 

(a) 


k 




? k 

\ 



FIG. 1: (Color online) (a) Energy spectrum of quasiparticles 
in graphene in the presence of circularly polarized irradia- 
tion (solid line) and in the absence of one (dashed line); (b) 
Scheme of quasiparticle transitions in graphene irradiated at 
the frequency cjq- 



energy gap ([7J between quasiparticle states originated 
from valence and conductivity bands of graphene. For- 
mally, the value of the photon-induced gap in the limit of 
large radiation intensities, given by Eq. ([7]), is the same as 
the value of dynamical gap induced at the Dirac point of 
graphene by a classical circularly polarized electromag- 
netic wave 10 ' 11 -. However, there is a conceptual differ- 
ence between the stationary energy gap ([7]) arisen from 
quantum-field theory and dynamical gapsiSr— arisen in 
graphene from classical fields. Historically, the term 'dy- 
namical gap' is come into the theory of graphene-field 
interaction-^— from the general theor y 14 i 15 developed 
to describe different quantum systems exposed to time- 
dependent classical fields and denotes the gap in spec- 
trum of Floquet quasienergies. Thus the dynamical gap 
is not true gap in the density of states of charge car- 
riers. As it was specially stressed in Ref. [lH, "since 
it is time-dependent problem, one cannot measure the 
gap directly from the density of states". In contrast to 
the dynamical gap, the gap ([7]) is obtained from station- 
ary solutions of the time-independent Dirac problem for 
the electron-photon system and is true gap in the den- 



sity of bound electron-photon states (charged quasipar- 
ticles). As a consequence, the stationary gap (UJ) will be 
to manifest itself directly in all phenomena sensitive to 
the density of states of charge carriers. In other words, 
the incorporation of such an additional physical factor as 
a quantum nature of electromagnetic field into the the- 
ory of graphene-field interaction, presented in the given 
paper, allows to consider the strong interacting electron- 
photon system in graphene as a gas of noninteracting 
quasiparticles with density of states defined by the en- 
ergy spectrum (|13[) . Since the quasiparticle concept is 
conventional approach of the modern physics to describe 
quantum systems with a strong interaction, the devel- 
oped theory is fruitful to predict and analyze different 
unexplored phenomena in graphene irradiated by light 
— particularly, novel physical effects discussed hereafter. 

If another electromagnetic field with the frequency w 
is applied to the considered electron-photon system, op- 
tical transitions of quasiparticles between its valence and 
conductivity bands can be possible. It follows from the 
energy conservation law that the transitions are allowed 
for uj > ujg, where uj g — e g /h. Therefore the gap (J7J re- 
sults in the threshold of absorption of an external field at 
the frequency u> — cu g . Since the threshold frequency, uj g , 
is controlled by the field intensity, Eq, this effect can be 
used as a basis for tunable detectors of electromagnetic 
radiation. The transition rate at the threshold frequency 
is defined by the intcrband dipole matrix element 

K^ (o)l^l^ (o))l = N^ 2 ;; ( + ^ o) i (i4) 

where d is dipole moment operator, and i = x,y. Besides 
optical transitions induced by the external field, there is 
the quasiparticle transition at the frequency luq, 

vfV o (k)^*+ o _ l(k ), (15) 

which is accompanied by one-photon absorption of the 
field ©. This transition, pictured in Fig. lb, needs a 
special discussion since the field © is inherent in the 
solved Dirac problem and cannot be considered as an ex- 
ternal field. It follows from Eq. (ITT]) that for u g < 
there is the wave vector, k', satisfying the condition 
£jY o (k') = e^ n _ 1 (k'). Therefore the transition (|T5|) is al- 
lowed by the energy conservation law for u> g < uiq and can 
take place at k = k' without the assistance of external 
energy sources. However the field ^ has been accounted 
exactly by the complete Hamiltonian of the electron- 
photon system ([3]) and cannot lead directly to transitions 
between different stationary states of the same Hamilto- 
nian. To realize the transition (fTSj) . the eigenstates of the 
Hamiltonian ©, ^)y n (k) and 5'^ o _ 1 (k), should be mixed 
by an external elastic scatterer (for example, by a defect 
of graphene lattice). Using Eq. (fT!?|) . we can calculate the 
matrix clement 

(*+ -r(k)|f/(r)|^ (k)), (16) 

which defines the rate of transition (fl"5"j) for any scat- 
tering field U(r). In the absence of scatterers, when 
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U(r) — (or, generally, U(r) = const), the matrix ele- 
ment (|16p is zero and the transition is impossible. Thus 
the quasiparticle transition (fl~5j) is scatterer-assisted pro- 
cess. 

It is seen in Fig. lb that excited quasiparticles, accom- 
plishing a cascade of non-radiative thermalization tran- 
sitions after the transition (|15[) . can create a population 
inversion between the conductivity band and the valence 
band at k = 0. Since the interband dipole moment ([14"]) 
is nonzero, that results in light emission from graphene 
at the threshold frequency, uj gi controlled by the field in- 
tensity, Eq. Therefore the scheme, pictured in Fig. lb, 
forms a physical basis for creating tunable sources of 
electromagnetic radiation. The proposed mechanisms for 
tunable detection and generation of electromagnetic ra- 
diation are especially challenging for the THz frequency 
range. Indeed, search for effective THz sources and de- 
tectors is one of most excited problems of modern ap- 
plied physic o 16 ' 17 . Moreover, the latest trend is the us- 
ing of different nanostructures to fill the THz gapi£~— . 
Therefore the application of graphene for THz devices fits 
well current tendencies in the graphene-based quantum 
electronics^ and the nanophotonics as a whole. 

Formally, the energy spectrum (|13p is of dielectric type 
for any field intensity, Eq. However the interband tran- 
sition (|15p leads to generation of free charge carriers in 
valence and conductivity bands for uj g < ujq. To turn 
graphene into true insulator, we need to forbid the transi- 
tion (|15[) by the energy conservation law. This forbidding 
corresponds to ui g > loq. Therefore the metal-insulator 
transition in graphene is a threshold effect and occurs at 
the critical field intensity, Eq = E c , satisfying the condi- 
tion Eg = hiOQ. Using Eq. |(7J), we obtain 



E r 



2v F \e\ 



(17) 



Being in the insulator state for Eq > E c , graphene has 
zero conductivity for the temperature T = and behaves 
as a semiconductor for T =/= 0. Applying the conventional 
theory of semiconductors 2 ^ to graphene for Eq > E c and 
T < e s , we can write the total density of free charge 
carriers for a single valley of graphene and for a certain 
direction of electron spin in the form 



m*T 



exp 



V 2T 



(18) 



where m* is effective mass of quasiparticles, which defines 
their energy spectrum, 



Ae(k) = ±e g /2 ± h 2 k 2 /2m 



(19) 



for small quasiparticle wave vectors k. Considering k- 
terms in the Hamiltonian as a perturbation and using 
the standard perturbation theory^, we can find the en- 
ergy spectrum of quasiparticles near their band edge as 
a k-power series expansion without assuming weakness 
of the field (the parameter of electron-field interaction, 



a = Wq/Hujq, is allowed to be not small). Obtaining 
terms ~ k 2 in the spectrum, we arrive at the exact ex- 
pression for the effective mass 



2e g {e g + 2frw Q )(s g + Fujq) 2 
v 2 F [{e g + 2hLOQf + e g ] 



(20) 



It should be noted that the energy spectrums of quasi- 
particles ([13]) and ([19]) are mutually complementary: 
Eq. (1131) is applicable within a broad range of wave vec- 
tors but only for small field intensities (a <C 1), and 
Eq. ([19]) with the effective mass ([20]) is valid for any field 
intensities but only for small wave vectors (HvFk/e g <C 
1). As expected, Eqs. (fT3"]) and ([19]) are equal for small 
field intensities and small wave vectors. 



IV. DISCUSSION AND CONCLUSIONS 

Finalizing the paper, we have to estimate the ef- 
fects discussed above. The distance between neighboring 
electron-photon states (j6]) is characterized by the two en- 
ergies, hujQ and e g = hu) g . Therefore the destructive in- 
fluence of scattering processes on quasiparticles can be 
neglected if uiqt » 1 and oj g T » 1, where r ~ 10~ 12 c 
is mean free time of charge carriers in graphene 2 -. Thus 
the developed concept of quasiparticles is applicable if 
the field ^ is both high-frequency and strong. Particu- 
larly, such a field can be created by using lasers. Let the 
wave is generated by a laser with the wavelength A and 
focused into a narrow beam with the diameter d ss A. 
Then the wave amplitude is Eq ps 4y/P/c\ 2 , where P is 
output power of the laser. Within this scheme the con- 
dition of metal-insulator transition (ITT]) can be satisfied, 
for instance, by using an ordinary low-power CCvlaser 
with A = 10.6 /im and P ~ 76 Watts. In this case the 
gap ([7]) is e g ~ 117 meV, that allows to observe the 
metal-insulator transition at room temperatures. As to 
the effect of THz emission from irradiated graphene, it 
takes place for substantially smaller powers P ~ 10~ 2 W. 

As it follows from the estimations given above, the con- 
sidered phenomena can be observable in graphene un- 
der a low-power irradiation. This is surprising, since a 
high-power laser radiation is usually necessary to observe 
effects of strong electron-photon couplin g 25 ' 26 . Let us 
give simple physical reasons to clarify this unexpected 
result. From the quasi-classical viewpoint, the charac- 
teristic energy of electron coupling to circularly polar- 
ized electromagnetic wave is the kinetic energy of elec- 
tron rotational motion induced by the wave. In the case 
of an usual condensed matter with parabolic electron en- 
ergy spectrum, e(fc) = h 2 k 2 /2m* , this kinetic energy is 
Wq — (cEq/uq) 2 /2m* . As to graphene, due to its lin- 
ear electron energy spectrum, s(k) — ±Hvp \k\, the en- 
ergy of the electron rotational motion induced by the 
wave is Wq = 2vfgEq/ujq. Since we have Wo/Wq ^> 1 
for weak fields Eq, the electron coupling to a low-power 
circularly polarized electromagnetic radiation is substan- 
tially stronger in graphene than in other condensed mat- 
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ters. The observability of the discussed effects for a low- 
power laser pumping is crucial from experimental view- 
point, since high-power lasers fiuidize a crystal lattice and 
are inapplicable to condensed-matter experiments. Thus 
graphene gives an unique opportunity to observe effects 
of strong electron-photon coupling, which cannot be ob- 
served in other condensed matters. As a result, the pre- 
sented theory opens a significant new area of graphene- 
related research, where condensed matter physics and 
quantum optics meet. 

We have demonstrated that the Dirac electron spec- 
trum of graphene leads to the strong electron interac- 
tion with circularly polarized photons for low light in- 
tensities. This results in bound electron-photon states 
(charged quasiparticles) which should be considered as a 



substantially new kind of field-matter coupling. The en- 
ergy spectrum of the quasiparticles is of dielectric type, 
that leads to the photon-induced metal-insulator transi- 
tion. From applied viewpoint, this effect can be used as 
a basis for creating tunable optoelectronic devices. 
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